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We examine the stability of matter-antimatter molecules by reducing the four-body problem into a 
simpler two-body problem with residual interactions. We find that matter-antimatter molecules with 
constituents (m±, m^~, m^~, fh±) possess bound states if their constituent mass ratio m\/m2 is greater 
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I. INTRODUCTION 

The study of the stability and the properties of matter-antimatter molecules has a long history, starting with the 
pioneering work of Wheeler who suggested in 1946 that (e + e~) might be bound with its antimatter partner to form 
a matter-antimatter molecule Since then, the problem has been examined theoretically by many workers [3-fl5| 
(for a review and other references see [IH). Wheeler went on to explore the properties of an assembly of (e + e - )™ 
atoms and molecules if they were made, and he outlined the phase boundaries in temperature and pressure separating 
various phases of (e + e _ ) atoms and (e + e~)-(e~e + ) molecules in their gaseous, liquid, super-fluid, crystal, and metallic 
states [17] . However, the experimental detection of matter-antimatter molecules is difficult, and the (e + e~)-(e _ e + ) 
molecule was successfully detected only recently, as late as 2007 [l||. 

In spite of extensive past investigations, our knowledge of matter-antimatter molecules remains rather incomplete, 
being limited to (e + ) m (e~) n and some aspects of (pe~)-(e + p). There are however stable and meta-stable charged 
particles and antiparticles, such as e^, pp, [i , w , K , t , etc. The conditions for the molecular binding of four- 
body particle-antiparticlc complexes containing these charged constituents are not known, nor are their annihilation 
lifetimes, if these matter-antimatter molecules turn out to be bound. 

With the advent of the Relativistic Heavy-Ion Collider at Brookhaven and the Large hadron Collider at CERN, 
a large number of charged particles and antiparticles are produced in high-energy pp and heavy-ion collisions. The 
production of matter and antimatter particles in close space-time proximity raises the interesting question whether 
chance encounters of some of the produced charged particles and antiparticles may lead to the formation of matter- 
antimatter molecules as debris of the collision. The detection of matter-antimatter molecules in high-energy nuclear 
collisions will need to overcome the difficulty of large combinatorial background noises that may be present. In another 
related area, recent production and trapping of cold antihydrogen |19h21| provide the possibility of bringing matter 
and antimatter atoms close together. Furthermore, electron-positron colliders with fine energy resolutions may be 
used to produce stable matter-antimatter molecules as resonances with finite widths, when the colliding e + and e~ 
combination has the same quantum number as the matter-antimatter molecules. 

The detection of new matter- antimatter molecules is however a difficult task, as evidenced by the long span of time 
between the proposal and the observation of the (e + e - )-(e~e + ) molecule. Additional instrumentation and experi- 
mental apparatus may be needed. It is nonetheless an interesting theoretical question to investigate systematically 
the general factors affecting the stability of matter-antimatter molecules, whether any of these matter-antimatter 
molecules may be bound, and if they are bound, what are their binding energies, annihilation lifetimes, and other 
characteristics. Answers to these questions will help us assess whether it may ever be feasible to detect them experi- 
mentally in the future. 

Following Wheeler [l[, we shall use the term "atom" to represent a two-body bound state of a positive and a 
negative charged pair that can form a building block, out of which more complex "molecules" can be constructed. In 
the present work, we shall limit our attention to molecules in which the four constituents {mi, mi, m%, 7714} consist 
of two charge conjugate pairs, with 7713 the charge conjugate of m-i, and 7114 the charge conjugate of mi. We shall 
arrange and order the constituents according to their masses such that mi>mi, and the charges of mi and mi be +e 
and — e respectively To make the problem simple, we shall consider molecules containing non-identical constituents 
mi and 7712 and their antiparticle counterparts, such that there are no identical particles among the constituents. 
Molecules with identical constituents require additional considerations on the symmetries of the wave function with 
respect to the exchange of the pair of identical particles, which are beyond the scope of the present investigation. 

To study the structure of the molecules, we shall consider only Coulomb interactions between particles and neglect 
strong interactions, as the range of strong interactions is considerably smaller than the Bohr radius of the relevant 
particle-antiparticle system. 

Previously, based on the method proposed for the study of molecular states in heavy quark mesons [22| , we obtained 
the interatomic potential for the HH system [23| . We shall generalize our consideration to cases of constituent particles 
of various masses and types and shall quantize the four-body Hamiltonian to obtain molecular cigenstatcs of the four- 
body system. If molecular states are found, we shall determine their annihilation lifetimes and their spin dependencies, 
if any. 

It is worth pointing out that the subject matter of molecular states appears not only in atomic and molecular 
physics, but also in nuclear physics and hadron spectroscopy. Wheeler's 1937 article entitled "Molecular Viewpoints 
in Nuclear Structure" introduced molecular physics concepts such as resonating groups and alpha particle groups to 
nuclear physics [24| . Indeed, nucleus-nucleus molecular states have been observed previously in the collision of light 
nuclei near the Coulomb barrier [25j . Molecular states of heavy-quark mesons have been proposed in high-energy 
hadron spectroscopy I22T [26Tj29| to explain the narrow 3872 MeV state discovered by the Belle Collaboration [3(| and 
other Collaborations [31[. The general stability condition established here for the Coulomb four-body problem for 
matter-antimatter molecules may have interesting implications or generalizations in other branches of physics. 

This paper is organized as follows. In Section II, we review the families of states of the four-body matter-antimatter 
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system so as to introduce the method of our investigation. In Section III, the mathematical details of our formulation 
are presented and the four-body problem is reduced to a simple two-body problem in terms of the interaction of 
two atoms with residual interactions. The interaction potential is found to consist of the sum of the direct potential 
Vdi r and the polarization potential V po i. In Section IV, we show how to evaluate the interaction matrix elements. 
In Section V, we show the analytical result for the direct potential Vdir- In Section VI, the polarization potential 
is evaluated for the virtual excitation to the complete set of bound and continuum atomic states. The results of 
the interaction potential for (pe + )-(e + p) are discussed in Section VII. The interaction potentials for other molecular 
systems are examined in Section VIII. In Section IX, we solve the Schrodinger equation for molecular motion and 
obtain the molecular cigenstates for different systems. We discuss the annihilation rates and lifetimes of the molecular 
states in Section X. Finally, Section XI gives some discussions and conclusions of the present work. Some of the details 
of the analytical results are presented in the Appendix. 

II. FAMILIES OF FOUR-PARTICLE STATES 

The four constituent particles can be arranged in different ways leading to different types of states. There is the 
AA family of states of the type j4(m^~m^)-^4(mjmj~), in which m\ and m^ orbit around each other to form the 
A„(m^~m^~) atom in the v state, while ffi^ and mj~ orbit around each other to form the A v i (mj m i ) atom in the v' 
state. When A v and A v < are separated at large distances, their asymptotic state energy is 

2 m\m,2 ( 1 1 \ 

E AA {n v ,n v >) = -a - — + — . (1) 

rri! + m 2 \nf, n„, J 

There is another MM family of states of the type M(m^fhi)-M(m2m2), in which mf and mj~ orbit around each 
other to form the M 1/ (mi~m^~) atom in the v state, while rfi^ and orbit around each other to form the M v i (m^ m^) 
atom in the v' state. Their asymptotic state energy is 

1 o / mi m,2 \ 

E MM\ n vi n v') = --, a — + • ( 2 ) 

For the same values of n v = n v i , the asymptotic state of the MM family lies lower in energy than the asymptotic 
state of the AA family, except when mi — 1112 for which they are at the same level, 

and E M M{n v ,n v ) = E A A{n v ,n v ) if m 1 =m 2 - (3) 

On the other hand, by varying the principal quantum numbers, many of the asymptotic states of one family can lie 
close to the energy levels of the asymptotic states of the other family. Level crossing between states and the mixing 
of states of different families can occur when the atoms are brought in close proximity to each other. 

There are two different methods to study molecular states. In the first method, one reduces the four-body problem 
into a simpler two-body problem. One breaks up the four-body Hamiltonian into the unperturbed Hamiltonians of 
two atoms, plus residual interactions and the kinetic energies of the atoms. The unperturbed Hamiltonians of the two 
atoms can be solved exactly. Using the atomic states as separable two-body basis, one constructs molecular states with 
the atoms as simple building blocks and quantize the four-particle Hamiltonian [221 ] . The quantized eigenstate obtained 
in such a method may not necessarily contain all the correlations. They may also not necessarily be the lowest states 
of the four-body system. They however have the advantage that the center-of-mass motion of the composite atoms 
are properly treated and the formulation can be applied to systems with vastly different mass ratios mi/m2- They 
provide a clear and simple picture of the molecular structure. They also provide vital information on the condition of 
molecular stability and the values of molecular eigenenergies. The knowledge of the molecular eigenfunctions provides 
information on other properties of the molecular states and their annihilation lifetimes. Furthermore, these molecular 
states can form doorways for states of greater complexity with additional correlations. For example, one can multiply 
the four-body wave functions of a molecular state obtained in such a method (see Eq. (fTQ|) below) by a complete set of 
correlated wave functions for a particular pair of constituents, and diagonalizc the four-body Hamiltonian with such 
a basis. The eigenstates obtained after such a diagonalization represent the splitting of the doorway state into finer 
molecular states containing additional correlations. As the method exhibits a clear molecular structure in terms of 
composite two-body objects, its asymptotic states illustrate how the molecule may be formed by the collision of the 
composite atoms. Finally, if one depicts the orbiting of one particle relative to another particle as a "dance pattern" 
with the topology of a ring, then the dance patterns of the different constituents and atoms in different families have 
distinctly different connectivities and topological structures. The transition of a state from one family to another 
family will involve the breaking of one type of dance pattern and re-establishing another type of dance pattern. It 
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is reasonable to conjecture that their distinct topological structures may suppress the transition amplitude between 
families and may allow the atoms to retain some of their characteristics and stability in their dynamical motion and 
transitions. 

There is an alternative second method to study the molecular states by taking the interaction potential to be the 
adiabatic potential obtained in a variational calculation for the lowest-energy state of the four particle complex in 
the Born-Oppenheimer potential, in which the positions of two heavy constituents are held fixed [8l-ll5j. It should be 
realized that such variational calculations have not yet been fully variational. By fixing the positions of the heavy 
constituents (the proton and antiproton in the case of the (pe~e + p) complex), the trial wave functions of the heavy 
constituents have been constrained to be delta functions without variations. If the motion of the heavy constituents 
were allowed to vary in a fully variational calculation, the lowest energy state would be the one in which the heavy 
constituents would orbit around each other in their atomic orbitals, with binding energies proportional to their heavy 
masses. The variational calculations of the higher molecular states will need to insure the orthogonality of the state 
relative to lower lying ones. Furthermore, as the variational calculation searches for the state with the lowest energy, 
motion in the relative coordinates between the heavier masses corresponds to constraining a trajectory along a path 
with adiabatic transitions whenever a level crossing occurs. In a collective molecular motion, these level crossing 
may not necessarily be adiabatic [HI because the speed of the collective motion becomes large at small interatomic 
separations, and the transition matrix element between different families may be suppressed due to the difference 
in their topological structures. It is reasonable to suggest that while variational calculations may provide useful 
information on the four-particle complex, they should not be the only method to examine the molecular structure of 
the four-body system. 

We shall use the first method to study the eigenstates of the four-body system and construct states built on the 
AA family. Accordingly, we break the four-body Hamiltonian into a part containing the unperturbed Hamiltonians 
of the A and A atoms and another part containing the residual interactions and the kinetic energies of A and A. 
The quantization of the four-body Hamiltonian then gives the molecular states of interest, as in a previous study of 
molecular states in heavy mesons in hadron physics [22J. Such a separation of the four-body Hamiltonian is justified 
because the composite A and A atoms are neutral objects, and their residual interaction Vj between the constituents 
of A and A involve many cancellations. As a consequence, the non-diagonal matrix element for excitation arising from 
the residual interaction is relatively small in comparison with the difference between unperturbed state energies, and 
the perturbation expansion is expected to converge. 

Molecular states can be constructed using composite objects of A v and A v i in various asymptotic v and v' states. 
We shall be interested in molecular states in which the building block atoms A and A are in their ground states 
at asymptotic separations. As the two ground state atoms approach each other, they will be excited and polarized 
and their virtual excitation will lead to an interaction potential between the atoms, from which the eigenstates of the 
molecule will be determined. 

Molecular states constructed with excited A v and A v i atoms can be similarly considered in a simple generalization 
in the future. Such a possibility brings into focus the richness of states in the four-body system, as the molecular 
states will have different interatomic interactions, obey different stability conditions, and have different properties 
with regard to annihilation and production. 



III. THE FOUR-BODY PROBLEM IN A SEPARABLE TWO-BODY BASIS 



In order to introduce relevant concepts and notations, we shall review the formulation of the four-body problem in 
terms of a simpler two-body problem |22j . We choose the four-body coordinate system as shown in Fig. [T] and label 
constituents m\ , 171% , to^~, and as particles 1, 2, 3, and 4, respectively with m,\>vni. The Hamiltonian for the 
four-particle system is 

j = l j = l k>j j=l 

in which particle j has a momentum p - and a rest mass rnj. The pairwise interaction Vjk(rjk) between particle j and 
particle k depends on the relative coordinate between them, rjh = rj — r&. 

We introduce the two-body total momentum Pj k = Pj +P k i and the two-body relative momentum p jfe = fkPj — fjPk j 
where 

„ m k 

Jk — , 

m 3 k 

rrijk = nrij + m k . (5) 
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A(34) 



FIG. 1. The coordinates of the {m^,m 2 , m^,m 4 } system. 



We choose to partition the four-body Hamiltonian in Eq. ((4]) into an unperturbed hi 2 + ^34 for atoms A(12) and 
A(34). a residual interaction Vr, and the kinetic energies of A and A, 

H = h 12 + h M + V I + 1 ^- + ^, (6) 
Imvx 2m 34 

where 



'»jfc = ^- + «jfc(ri-r fc ) + m jfc for {jfc} = {12} and {34}, (7) 
and 

Vr = wi4(ru) + vi 3 (r 13 ) + v 23 (r 23 ) + v 2 4,(r 2i ). (8) 

The eigenvalues of the unperturbed two-body Hamiltonians hi 2 and h 34 can be solved separately to obtain the bound 
state wave functions and masses Mj).(v) of atoms A and A, 

h jk \(jk) v ) = [e jk {v) + m jk ]\(jk) v ) = M jk {v)\{jk) v ). (9) 
When we include Vj as a perturbation, the eigenfunction of H becomes [33| 

* (r,r 12 ,r 34 ) = ip(r)< \A A ) - > — — — (10) 

where r = R\ 2 — H34 is the interatomic separation (see Fig. 1), Rj k = fj rj + f k r k is the center-of-mass coordinate 
of mj and m k , and indicates that the sum is over a complete set of atomic states |A^Av), including both bound 

and continuum states, except |^4o^o)- The eigenvalue equation for the four-body system with eigenenergy e is 

HV(r,r 12 ,r 34 ) = [M 12 (0) + M 34 (0) + e]tf(r, n a , r 34 ). (11) 

Working in the center-of-mass frame and taking the scalar product of the above equation with |j4oAo), the eigenvalue 
equation for the four-body system becomes the Schrodinger equation for the motion of A (12) relative to A (34), 

£- +V(r))i>(r) = e^(r), (12) 

where p is the relative momentum of the composite particles (atoms) 

_ M 34 (0)Pi 2 -M 12 (0)P 34 ml 
P Mi 2 (0) + M 34 (0) ' [ ' 



fJ-AA is the reduced mass of the two atoms 



_ M 12 (0)M 34 (0) h& , 
~ Mi 2 (0) + M 34 (0) ' U4j 



and the interaction potential V(r) in Eq. (fT2|) is given by 



IAA\V\AA\ S^' \{A^\Vi\MA^ 
V{r) = {AqAo \Vi \A Ao) - > ttt- 7T77 77- — -. (15 

'A,A' e j4 (A) + e^A') - ca(0) - e^(0) 
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We call the first leading-order term on the right hand side of the above equation the direct potential, Vdi r (r), 

V dil (r) = (AoAoIV^AoAq), (16) 

which represents the Coulomb interaction between the constituents of one atom and constituents of the other atom. 
We call the second next-to-leading order term in Eq. (fT5|) the polarization potential, V po i(r), 

V ^ r > l^yy eA (X) + e A {\>) e A (0) - e A (0) ' U '> 

which is always negative. It represents the effective interatomic interaction arising from the virtual Coulomb excitation 
of the atoms as they approach each other. 

To study molecular states based on A and A atoms as building blocks, we quantize the Hamiltonian for the four- 
body system by solving the Schrodinger equation (|12[) . For such a purpose, our first task is to evaluate the interaction 
potential V(r) in ([T5"j) by calculating the direct and polarization potentials in (fT6"| and ([TT)) . 



TABLE I. The atomic length unit, Bohr radius, a — h/afi, and the atomic energy unit, a 2 /x, for different combinations of 
atomic constituents. 









K+ 


P 


r+ 


e 


a (fm) 
a 2 n (eV) 


53174 


53111 


52972 


52946 


52932 




27.08 


27.11 


27.18 


27.20 


27.20 




a (fm) 




449.7 


310.7 


284.7 


271.1 




a 2 /i (keV) 




3.202 


4.635 


5.057 


5.311 




a (fm) 






248.5 


222.6 


209.0 




q 2 /i (keV) 






5.794 


6.47 


6.891 


K~ 


a (fm) 
a 2 /i (keV) 








83.59 
17.22 


69.99 
20.57 


P~ 


a (fm) 
a 2 /i (keV) 










44.04 
32.69 



For molecular states in the AA family, the atomic units of A and A are the same. To exhibit our results, it is 
convenient to use the atomic unit of the ^(m^m^") (or A) atom as our units of measurement: 

1. All lengths are measured in units of the Bohr radius of the ylfro^m^) system, ai 2 = ^0^12, where [i\2 = 
mim 2 /(mi + m 2 ). 

2. All energy are measured in units of a 2 /xi 2 for the A(m^ m^) system, which is two times the Rydberg energy, 

2£Ryd{12}- 

3. As a consequence, the reduced mass [i AA in the Schrodinger equation (fT2|) for molecular motion in coordinate 
r needs to be measured in units of /112, 



{n AA in atomic units) 



(J-AA _ A/ 12 (0)A/ 34 (0) _ (rm +m 2 ) 2 
M12 Ati2[Mi 2 (0)+M 3 4(0)] ~ 2m 1 m 2 



(18) 



For brevity of notation, we shall omit the subscript {12} in ai 2 and a 2 /zi 2 except when it may be needed to resolve 
ambiguities. In Table I, we show the physical values of the Bohr radius a and the energy unit cv 2 /z. They are different 
for different ^4(and ^4) atoms that build up the ^4^4 molecule. These quantities will be needed in Section X to convert 
atomic units to physical units. 



IV. METHOD TO EVALUATE (A x A y \v jk (r jk )\A A ) 

To obtain the direct and polarization potentials, we need to evaluate the matrix element (A\Ay\Vjk(rjk)\AoAo). 
We shall examine first the case when both A\ and Ay are bound states. The case when one or both of A\ or Ay 
lies in the continuum necessitates a different method and will be discussed in Section VLB. 

When both A\ and Ay are bound, we can use the Fourier transform method to evaluate the matrix elements 
dHHi)]- Here, A\(12) and Ay (34) can be represented by normalized hydrogen wave functions (j> x {r 12) and ^(^34), 



respectively. The residual interaction Vj k (rj k ) is a function of rj k . We need to express Vj k in terms of r, r 12 , and 

^34, 

r jk = Tj - r k =r + F A (jk) r 12 + F A (jk) r 34 , (19) 

where the coefficients F^ AA y(jk) have been given, with a slight change of notations, in Ref. (22^ . 

F A (U)= / a , F A (U)= f 3 , 

F A (13)= f 2 , F A (13) = -.U, 

F A (23) = -fi, F A (23) = -/ 4 , 

F A (24) = -/i, F*(24) = / 3 . 

The matrix element (A^Av li^rj/^jAoAo) can be written as 

(A x Ax>\v jk (r jk )\AoAo) = J dr 12 dr 34 pxo(r 12) Pyo( r 34)v ]k (r + F A (jk)r 12 + F A (jk)r 34 ), (20) 
where p A {ri2)=4>* x (ri2)4>o( r i2) and Pvo( r 3 4 ) =< / > A' ( r 34)<fo( T "34)- Introducing the Fourier transform 

P A xo A (P) = J dyejP y p A XQ A (y), (21) 

and 

we obtain 

(A x A x ,\v jk (r jk )\A A ) = f J^ e ^ r p A (F A (jk)p)p A (i 



; (p) = J dr jk e-^ r ^v jk (r ]k ), (22) 



For our Coulomb potential 



«* . Sa, (25) 
where &j is the charge of rrij , the Fourier transform of the Coulomb potential is 

v jk(p) = —jp ■ ( 26 ) 

The Fourier transform p A (F A (jk)p) depends on the sign of F A (jk) and the I quantum number of \A\) = \nlm). It 
is easy to show that 

Pi {F A {jk)p) = [sign(F A (jk))} 1 p A (f A (jk)p), (27) 
where sign(F A (j k)) is the sign of F A (jk), and f A (jk) is the magnitude of F A (jk), 

f A (jk) = \F A (jk)\. (28) 

Substituting Eq. ([27]) in Eq. (|23]t. we obtain 



(A x Ay\v jk (r jk )\A A ) = s / ^eP- r &(f A Vk)p)tf, (f A (jk)p)v jh (pl (29) 



where s is the sign factor 



(27T) 

= [BignF A 0'*)]' [signF^O-fc)]'. (30) 
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FIG. 2. The direct potential Vdi T (r) in atomic units, for different systems with different mi/rri2. 



THE DIRECT POTENTIAL V di Jr) 



The direct potential is equal to the matrix element (j4 j4o| Vj |j4 j4 ) . We can apply the results of Eq. (|29[) to 
evaluate this matrix clement. We note that 



Poo(p) 



16 



((pa) 2 +4) 2 - 



Substituting this into Eq. (j2U)) . we have 
(A A \v 3k (r jk )\A A ) = 



dp 



oipr. 



16 



16 



47ra 



((p/^) 2 +4)2((p/^a 3 4) 2 +4)2 p 2 



(31) 



(32) 



which leads to the direct potential 



V dix (r) 



ae 



-2r/f ia 



r 

1 r 
'4M 



+ - 



ae 



-2r// 2 a 



1 r 
4/ 2 a 



VtUt ~ 3/|) 

(A 2 -/!) 3 

r \ 1 r 
Ja~) + 24 M 

2/K/l - 3/ 2 
(/ 2 2 -/i 2 ) 3 
1 r 



2\2 



r 
A a 



(A 2 - /I) 



6— + 3 



r 
A a 



1 



2/| 



ha) 24 / 2 a 



r 
/2 a 



(/ 2 2 -/i 2 ) 2 
r 

6— + 3 

ha 



r 
ha 



(33) 



This direct potential is a sum of two Yukawa potentials of screening lengths /ia/2 and / 2 a/2, multiplied by third-order 
polynomials in r. It can be shown numerically or analytically that for mi = m 2 , the direct potential Vd; r is zero. 
In the region close to r — >■ 0, the direct potential becomes 



lim Vdir(r) 

r->0 



(A 2 -/!) 3 

ae -2,// 2 a / 2/ 4 (/ 2 _ 3/ 2) 

(A 2 -A 2 ) 3 



-2r//ia 



A a 

ae -2r// 2 a 



2A 4 



2\2 



A a 



(A - A) 
f 2 4 



2/ 2 4 



(A 2 -A 2 ) 



2^2 



If mi 3> m 2 in this region close to r — > 0, then the direct potential becomes 



lim Vdir(r) 

r— mi>m2 



ae 



-2r/f ia 



ae 



-2r/f 2 a 



3a 



■ r )n 



4/i a 4/ 2 a 



(34) 



(35) 



9 



For this case of mi 3> m 2 , we have /ia/2~a mi7712 /2 and /2a/2~a mi m 1 /4, where a mjTO .=(mj+mj) /amjTOj is the 
Bohr radius of the (mirrij) system. The first term is a screened Coulomb interaction with the range of a miril2 /2 and 
the second term is a repulsive screened Coulomb potential with a range of o mi#il /4. The last three terms reduce the 
binding energies of molecular bound states in the case of mi 3> ui2- 

Equation Q33p is a general result applicable to any mass ratio of mi/m^, and is a generalization of the results of [3[ 
that represents only an approximation for mi/m^ 3> 1. 

We show in Fig. 2 the direct potential in atomic energy units, a 2 /i, for (pe~)-(e + p), (p/j~)-(/j + p), and (r + p)-(pr _ ), 
as a function of the interatomic separation r in atomic units, a. For systems with a large ratio of mi/ mi, the 
interaction is slightly repulsive at large separations, owing to the repulsion of like charges. The repulsive interaction 
of like charges is strongest when the atoms are nearly "touching" each other at r ~ 2a, leading to development of 
a barrier there. At r < a, the interaction between the heavy unlike charges dominates, and the direct potential 
changes to become strongly attractive. We observe in Fig. 2 that as the mass ratio mi/m2 approaches unity, there 
is a cancellation of both the attractive and repulsive components of the direct potential. The repulsive barrier is 
lowered and the direct potential becomes less attractive at short distances. In fact, as we noted earlier, Vdir vanishes 
if mi = I7i2. 



VI. THE POLARIZATION POTENTIAL V pol (r) 

The polarization potential Vp Q i(r) is the effective interaction between A and A arising from virtual atomic ex- 
citations. It can be obtained as a double summation over A\ and Ay involving the excitation matrix element 
(A\Ay\Vjk(rjk)\A A ). The summation over A\ and Ay states includes the complete set of bound and continuum 
states, but excludes the ground states. We shall make the assumption that the virtual excitation is predominantly 
electric dipole in nature and shall truncate the set of excited A\ and Ay states to include only I = 1 states. Because 
the ground states Aq and Aq have no orbital angular momentum, the azimuthal quantum numbers m of A\ and 
Ay must be equal and opposite. The polarization excitation therefore contains contributions where A- A' are bound- 
bound (bb), bound-continuum (be), and continuum-continuum (cc), with 1=1. We shall discuss separately how these 
excitation matrix elements can be evaluated. 



A. Bound-bound excitation matrix elements 



When both A\ and Ay are bound states, the "bound-bound" excitation matrix element (A\Ay \Vi\AqBq) can be 
evaluated using the method of Fourier transform. The results were presented previously in [23| ■ We shall rewrite the 
same result in a slightly simplified form. As shown in Appendix A, the relevant matrix element for \A\) = \nlm) and 
\A\t) = \n'l {—m)) with I = 1 is given by 

(A x By\v jk (r jk )\A B ) = s J ^R* lm>W0 [fAWp] R£, H - m)tW0 \fAWp] 

x Vjk(p)J(p,r), (36) 
where s is the sign factor as given by (f3T))) . J{p,r) is given in terms of the spherical Bessel functions, 

T( s j jo{pr) - 2j 2 {pr) for m = 0, 

{-[jo{pr) + J2[pr)\ for m = 1, 

-Rnim.iooO?) is given in terms of the Genegbauer polynomial C", 



WW - m 3 £ ( n ~l + 3 ) /r— (i _ 0) 



n+l V2/<H,V K 



K=Q 



npa(n+l) 2 2 6 (n- K ) 2 ( {npaf - (n + l) 2 

2 \3 Ly ri-2-K \ i ^2 , f m , 112 J ' 



{{npa) 2 + (n + l) 2 ) 3 n ~ A ~ K \ (npa) 2 + (n + 1) 

the variable 8 is 



(39) 
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and the normalization constant N n i is 



y a J ri 4 [(n + Z)!J 

Thus, the six dimensional integral over r\2 and in the matrix element is reduced into a one-dimensional integral 
that can be readily carried out numerically. 

B. Bound-continuum and continuum-continuum excitation matrix elements 

For a given interatomic separation r, the excitation matrix element (A\ A\i \ Vi\AqBq) involving one or two continuum 
states can be evaluated by direct numerical integration in the six-dimensional space of r 12 and r 34 . As constrained by 
the ground state wave functions of Aq and Aq, the integrand in such an integration has weights concentrated around 
the region of ri2 ~ and r 3 4 ~ 0, and the continuum wave function does not need to extend to very large distances. 

Following Bcthe and Salpeter [35|], we use the radial wave function for a continuum state of A (or A) with momentum 
\k\ = k as given by 

R k i(r) = ^x^F^kr), (41) 
kr V 7r 

where Fi(r),kr) is the regular Coulomb wave function (jtj, Hj. The coefficient of the wave function has been chosen 
according to the normalization 



i in 



k 2 dk 



Y r f-F^ kr)Y lm {6, 0) \ ( 1 Jl Fl („, kr)Y lm (9, 



(42) 



For the excitation to a continuum state in A and a bound state in A this closure relation allows us to write the 
bound-continuum contribution to the polarization potential to be 

vbc(r) _ y /' KRiir^YUOMAylVjlAoAo^ 

There is a similar contribution for the excitation into a bound state in A and a continuum state in A. 
The continuum-continuum contribution to the polarization Vp^(r) is given similarly by 

vu r ) = -J2 J k2dk J2 J k ' 2dk ' 

lm V 

K^(ri2)^ m ^12</'l2)i?^,(r34)^- m (034034)|^|AoA o )| 2 



e A (k) + e A (k') - e A (0) - e A (0) 



(44) 



To evaluate the excitation matrix element of Vi in Eqs. (|43|) and (|4"4")) . we discretized the continuum momentum k 
(and k') into momentum bins. For each of the bins, the wave functions in terms of the the spatial coordinates r\i and 
r34 are all known. We shall again limit our consideration to dipole excitations with I = 1' = 1 only. In the numerical 
calculations, the residual interaction Vj is a sum of four Coulomb interactions which depend on the magnitude of the 
radius vector where {jk} = {14}, {13}, {23}, and {24}. The square of the magnitude \rjk\ 2 can be evaluated as 

\r jk \ 2 = r 2 +r\ 2 + r| 4 + 2F A (jk)rr 12 cos Cl(r,r 12 ) 

+2F B (jk)rr 3i cos n(r,r 3i ) +2F A (jk)F B (jk)r 12 r 34: cos Cl(r 12 ,r 3i )- (45) 

It is convenient to choose the coordinate systems of r, r\2, and r34 such that their z-axes lie in the same direction, 
and their corresponding x- and y-axes are parallel to each other. With this choice of the axes, we have 

cos Q(r, r i2 ) = cos 8±2, (46) 
cos (r , r 34 ) = cos 6> 34 , (47) 
cosf2(r 12 ,r 34 ) = cos0 12 cos6» 34 + sin6» 12 sin6» 34 cos(0 12 - (f> 34 ). (48) 

These relations allow us to determine the integrand and carry out the 6-dimcnsional integration in r 12 and r 34 , for 
the evaluation of the excitation matrix element and the polarization potential. 
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FIG. 3. The quantity r G V po i(r) for different polarization potential components plotted as a function of r in atomic units, for 
the (pe~)-(e + p) system. The numbers on the right are the Ce coefficients from Eisenshitz and London for H2 (38( [. 



VII. THE INTERACTION POTENTIAL IN {pe~)-(e+p) 



Using the methods discussed in the last section, excitation matrix elements can be evaluated and the polarization 
potential V po i from different contributions can be obtained. We show in Fig. 3 the quantity r 6 V po i(r) for the (pe~)- 
(e + p) system as a function of r, where r V po \(r) from the bound-bound, bound-continuum, and continuum-continuum 
contributions are shown as the dashed, dash-dot and dash-dot-dot curves respectively. In these calculations, we 
include bound states up to n — 20 in bound-bound calculations and n = 12 in bound-continuum calculations. For 
calculations with continuum states, we include states up to k = 6 atomic units. 

We note in Fig. 3 that the curves of r 6 Vjp i flatten out at large values of r. This indicates that the attractive 
polarization potentials behave asymptotically as — Cg/r 6 , the well-known van der Waals interaction at large distances 
between atoms. The asymptotic values of the Cq from different contributions have been give n along with the theoretical 
curves. They have also been obtained previously by Eisenschitz and London for H2 [38[ as listed on the right side 
of the figure. There is reasonable agreement between the Cg values obtained in the present calculation and those 
of (38|. There is a small difference between the Cq numbers involving continuum states with those of (38|. These 
small differences may arise from the fact that the curves involving continuum states have not yet become completely 
flattened and thus they may have not yet reached their asymptotic values. 

From these results, we note that at large separations, the bound-continuum contribution is much larger than the 
continuum-continuum contribution and is slightly smaller than the bound-bound contribution. 

The situation is different at small separations. We plot V po \(r) as a function of r in Fig. 4 for (pe~)-(e + p) and 
r < 5a pe . We find that the bound-bound contributions are smaller than the bound-continuum contributions, which in 
turn are smaller than the continuum-continuum contributions. The total polarization potential remains well-behaved 
at small r. Its magnitude is much smaller than the magnitude of the direct potential dominated by the attractive 
screened potential between p and p. 

Having obtained both the direct and the total polarization potential, we can add them together to obtain the 
interaction potential V(r) = Vdh(r) + V po \(r). We show in Fig. 5 the interaction potential V(r) and its components 
Vdir(r) and V P oi(r) for the (pe~)-(e + p) system. We note that for this case of large ratio of constituent masses m\/rri2, 
the polarization potential is small compared to the direct potential at short distances and the total interaction is 
attractive at r < a. The repulsive barrier at r ~ 2a that comes from the direct potential remains. The repulsive 
interaction decreases at larger separations. There is a pocket structure at r ~ 6a pe that is very shallow and arises from 
the interplay between the repulsion of like charges at intermediate distances and the attractive polarization potential 
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VIII. THE INTERACTION POTENTIAL IN DIFFERENT SYSTEMS 



It is of interest to see how the interaction potential and its various components vary as a function of the constituent 
masses. The various potential components and the total polarization potential for (/i + e~)-(e + /i~) in atomic units are 
very similar to those of (pe~)-(e + p) and will not be presented. The situation changes slightly for (pfi~)-(fi + p). We 
show the polarization potential components in Fig. 6 and the total interaction potential in Fig. 7 for (p^,~)-(^, + p). One 
notes from Fig. 6 that for this case with m p /m M = 8.88, the bound-bound contributions to the polarization potential 
dominate over the bound-continuum or the continuum-continuum contributions. The results in Fig. 7 indicate that 
the total polarization potential is however small compared to the direct potential at r < a. The other features of the 
interaction potential is similar to those of the (pe~)-(e + p) system. 

When the ratio of constituent masses mx/m^ approaches unity, the qualitative features of the different components 
change significantly. In Fig. 8, we show various components of the interaction potential for the {t + p)-(pt~) system, 
for which m T /m p = 1.89. As one observes, the polarization potential is dominated by the bound-bound component 
while the bound-continuum and continuum-continuum contributions are small. The direct potential is much reduced 
compared to the case with large constituent mass ratios m\/m2 and is now of the same order of magnitude as the 
polarization potential. 
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FIG. 7. The direct potential, polarization potential, and the total interaction potential for the (p/j, )-(^ + p) system. 



IX. QUANTIZATION OF THE FOUR-BODY SYSTEM 

To study molecular states based on A and A atoms as building blocks, we quantize the Hamiltonian for the four- 
body system by solving the Schrodinger equation (fT2]l. The states of the system depend not only on the interaction 
potential V(r) but also on the reduced mass. 

Because we use the atomic units of the A{m\m2) atom to measure our quantities, it is necessary to measure the 
reduced mass for molecular motion in units of /Z12, as given by Eq. (|18|) . 

/ • , • • j \ (™i +m 2 ) 2 (1 + mi/m 2 ) 2 

{Haa m atomic units) = — = — ; . (49) 

277117712 2mi/m2 

To provide a definite description of the constituents, we order the masses of the constituents such that m\>7U2 and 
characterize the system by the ratio m\jmi. 

In our problem, the use of the atomic units of A and A as described in Section III brings us significant simplicity. 
We have just seen that the reduced mass in atomic units is a simple function of the constituent mass ratio mi/7712. 
It should also be realized that the interaction potential and its different components in atomic units depend only 
on the various coefficients F AA or /1 2 3,4) which are themselves ratios and are uniquely characterized by mi/m,2- 
Therefore, the Coulomb four-body system in atomic units is completely characterized by mi/m 2 - Consequently, two 
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FIG. 8. The direct potential, polarization potential, and the total potential for the (r + p)-(pr ) system. 
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FIG. 9. The eigenenergies of the molecular states states in (pe )-(e + p), (K + e )-(e + K ),(ir + e )-(e + Tv ), and (fi + e )-(e + (i ), 
in atomic units. 



different four-body systems with the same mxjmi will have the same molecular state eigenenergies and eigenfunctions 
in atomic units. We can infer the stability of a molecule by studying the change of the molecular eigenenergies as a 
function of m\/m2- 

We give in Table II the values of mi/ni2, the reduced mass in atomic units of the A atom, and the molecular state 
binding property, for many four-particle systems. The reduced mass decreases as 771 1/771 2 decreases. 

How docs the interaction potential varies as 7711/7712 decreases? For the case of 7711/7712 3> 1, the interaction potential 
at small r is dominated by the direct potential over the polarization potential. The total interaction potential is 
strongly attractive at small r. As 7711/7712 decreases and approaches 1, the direct potential is significantly reduced, 
and the interaction potential becomes dominated by the polarization potential. The net result is a decrease in the 
strength of the attractive interaction. 

We solve the Schrodingcr equation (TT^I) to obtain the eigenstate energies e for different molecular systems using 
the corresponding reduced masses and interaction potentials. The fourth column in Table II indicates whether bound 
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FIG. 10. The eigenenergy of the Is molecular state as a function of roi/m2. 



TABLE II. Relationship between mi/7712, reduced mass, and the presence of bound molecular states in many four-body 
systems. 



System 




Reduced mass 


Eigenenergy of 


A(m|m2 )-i(m^"mj") 


mi/mi 


in atomic units 


Is molecular state 


(mi > 7T12) 








(pe~)-(e + p) 


1836.2 


919.1 


-306.2 


(K+e-)-(e + K-) 


966.1 


484.1 


-160.6 


e~)-(e + -k~) 


273.1 


137.6 


-44.63 


Gu+e-)-(e+/0 


206.8 


104.4 


-33.52 


(pfi )-(fi + p) 


8.88 


5.50 


-0.540 




4.67 


3.44 


-0.0239 


(t + p)-(pt ) 


1.89 


2.21 


None 



states are present in various systems. There are bound states in (pe~)-(e + p), (K + e~)-(e + K~), (7r + e~)-(e + 7r~), 
(/i + e~)-(e + yu~), (pu~)-(u + p), and (K + u~)-(u + K~). The results in Table II indicate that bound molecular states 
exist for four-particle systems if m\jmi is greater than about 4 (or if the reduced mass is greater than or about 3 
atomic units). The eigenenergies for many four- particle systems are shown in Fig. 9. Wc label eigenstates by the 
angular momentum and the principal quantum number, which is equal to the number of nodes plus I + 1. The 
eigenenergies (measured in their corresponding atomic unit a 2 /X12) move up into the continuum as mi/m2 approaches 
unity and the reduced mass decreases. 

We can examine the molecular state energies of (pe~)-(e + p) and compare them with those of a pp atom. For the 
(pe~)-(e + p) molecule, the molecular Is state energy is located at -306.2a 2 /i pe and the ns state is higher than the np 
state, whereas the atomic Is state energy for a pp atom is -459.5a 2 ^ pe and the ns state has the same eigenenergy as 
the np state. The differences in state energies and ordering arise because in the (pe~)-(e + p) system, the interaction 
potential at small r is given by Eq. (|3"3")) (or Eq. ([3"5]) ') that contains the first term of an attractive screened potential 
with a screening length f\a/2 ~ a/2. This screening length is so large compared to the pp Bohr radius that the 
attractive screened potential is nearly —a/r in character, as in a pp atom. However, there is an additional repulsive 
screened Coulomb interaction in Eq. (f33|) with a screening length /2a / 2 that is comparable to the pp Bohr radius. 
This repulsive screened Coulomb interaction arises because in the (pe~)-(e + p) molecule, the orbiting of the leptons 
with respect to the baryons leads to the motion of the baryons. As a consequence, the proton and the antiproton 
have a spatial distribution and the interaction between the charge distributions of the proton and the antiproton is 
reduced from their point-charged values. This additional repulsive screened Coulomb interaction in Eq. (|33|) raises 
the eigenenergy of (pe~)-(e + p) in the Is state molecular state relative to the eigenenergy of the pp atom in the Is 
state, and the ns state to lie higher than the np state. 

For (pix~)-(p + p) for which mi/7712 = 8.88, Is state energy is at -0.5425 atomic units. The molecular state is weakly 
bound. For (K + p~)-(p + K~) for which TOi/m 2 = 4.67, the Is state energy lies at -0.0239 atomic units , which is 
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just barely bound, indicating that mi/m2 ~ 4 is the boundary between the region of bound and unbound molecular 
states. For (r + p)-(pr _ ) for which m\/m 2 = 1.89, there is no bound molecular state. 

We plot in Fig. [10] the eigenenergy of the molecular Is state as a function of mi/m2. As one observes in Fig. [10] 
and Table II, the engenenergy varies approximately linearly as a function of mi/m2, with an intercept of E = at 
7711/7712 ~ 4, indicating that bound molecular states are present for mi/ mi greater than about 4. 



X. ANNIHILATION LIFETIMES OF MOLECULAR STATES 



Having located the energies of matter-antimatter molecules in various systems, we would like to calculate their 
annihilation lifetimes. The annihilation probability depends on many factors: the spatial factor of contact probability 
that is determined by the particle wave functions, the type of annihilating particles whether they are leptons or 
hadrons, and the spins of the annihilating pair if they are leptons. We shall discuss these different factors in turn. 



A. Spatial Factor in Particle-Antiparticle Annihilation 



In our matter-antimatter molecules, mi-7774 and 7712-7773 are charge conjugate pairs which can annihilate. The total 
annihilation probability naturally comprises of Pi 4 for the annihilation of mi and 7714, and P23 for the annihilation of 
m2 and 7713 . These spatial factors Pjk can be obtained by approximating the constituent wave function to contain only 
the first term of the perturbative expansion in Eq. (fit)]) , as amplitudes of the excited states relative to the unperturbed 
states are small and the excited states have greater spatial extensions that suppress the annihilation probabilities. 

We shall limit our attention to the annihilation of s-wave molecular states, which dominates the annihilation process. 
By the term "annihilation" in the present work, we shall refer to the annihilation in the s-wave molecular states only. 
The probabilities for the annihilation in higher angular momentum states are higher order in a [39| and involve not 
only derivatives of the molecular wave function at the origin but also complicated angular momentum couplings. They 
will need to be postponed to a later date. 

We shall calculate the spatial factor Pjk for the molecular (s-wave) state ip(r) built on A v (mim 2 ) and A;,/ (77737774) 
atoms. It is quantitatively defined as the probability per unit volume for finding the conjugate nij-mk pair to be 
at the same spatial location, in the molecular (s-wave) state ip( r ) with A v {m\m 2 ) and (77737774) atoms. It is the 
expectation value of 5(rj k ), 



P jk = J drdr 12 dr 3i At(r 12 )At,{r M )ij*(r)5{r jk )iP{r)A„(r 12 )A„,{r 3i ). (50) 
We shall consider molecular states built on v = v' = 0, then 

P ih = I dr\iP(r)\ 2 (A A \6(r jk )\A A ) 

dr\^{r)\ 2 (A A \S(r + F A (jk)r 12 + F B (jk)r M )\A A Q ) , (51) 



in which (AoAo\5(r jk)\AoAo) has the same structure as Vdh- except that Vjk(rjk) is replaced by a delta function. For 
this case with v = v' = 0, the sign of F A A does not matter, and we can just use f A A for F A A . Similar to Eq. (|32l) . 
we have 

(A Ao\S(r jk )\A A ) = [ ^e^ ^ - 1 f„ ■ — 2 ,, 1 f, ■ /| .„ . (52) 



(2tt)3 ((p/Aa ) 2 +4)2((p/^a ) 2 +4)2 

For both {jk} = {14} and {23}, f A (jk) - f A (jk), and we get 

D jk (r) = (A A \5 jk (r jk )\A A ) 
_ e~ XA ^ k > [X A {jk)f 



4tt 48 



{{X A (jk)r] 2 +3X A (jk)r + 3}, (53) 



where 



X A (jk) = 



2 



fA{jk)a 



/a(14) = /x(23) = 

mi + m 2 771l + 7712 
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Eq. (f5"Tj) becomes 

P jk = J dr\^(r)\ 2 D jk (r). (54) 

As the wave function ip(r) has been obtained from our solution of the Schrodinger equation, the spatial factors can 
be calculated numerically. 

Following Wheeler [l|, we can consider an antiparticle fh to be at rest while its conjugate particle m sweeps by 
with a annihilation cross section c™™ at a velocity v, clearing a volume er™™i> per unit time. The probability of 
finding the pair of particles j and k to be in contact, per unit spatial volume, is Pjk- As a consequence, the number 
of particle-antiparticle contacts per unit time, which is equal to the rate of annihilation wj^ m , is given by 

< = Ox^- (55) 



B. Annihilation of Lepton Constituents 



The matter-antimatter molecules we have been considering consist of both leptons and hadrons. The annihilation 
cross section cr™™ depends on the particle type and the total spin of the annihilating pair. We shall discuss the 
annihilation of lepton pairs in this subsection. The annihilation of hadron pairs will be discussed in the next subsection. 

The mechanism for the annihilation of a lepton pair is well known 0, |4(| . It proceeds through the electromagnetic 
interaction with the fine-structure coupling constant a and the emission of two or three photons. A lepton pair with 
spin 5=0 can annihilate only into two photons, and a lepton pair with spin S=l can annihilate only into three photons 
[l[ . The lepton pair annihilation cross section multiplied by velocity v is given by [l], |4^| 



a s jk =o v (f° r annihilation into 2 photons) = it ^ — ^ J (56) 

4(tt 2 - 9) / a \ 2 



\ 2 

a 



°s k =i v (f° r annihilation into 3 photons) = — - — - — — a ( — 2 ) (57) 

y y iic j 

As a consequence, the rates of annihilation of a lepton pair in the singlet and triplet states are P, H(J 



2 

a 



w S jk =o (annihilation into 2 photons) = tt I — - J Pjk (58) 

n , s 4 <> 2 -9) / «\ 2 

ws fe =i (annihilation into 3 photons) = a I — 2 I Pjk (59) 

y \ fie J 

Thus, to obtain the rate of annihilation of a lepton-antilepton pair {jk} in a molecular state, it is necessary to find 
out the probabilities for the pair to be in different Sjk pair spin states in the molecule. 

Our molecular states have been constructed by building them with A(12) and .4(34) atoms in their ground states. 
As we restrict our considerations to only s-wave molecular states, there is no orbital angular momentum between the 
A and A atoms. In picking the lepton I from atom A and the antilepton I from the other atom A, the probabilities for 
different lepton- antilepton pair spin states depend on the angular momentum coupling of the lepton-antilepton pair 
with the remaining constituents. 

We consider first the case when the remaining constituents are also fcrmions. The four constituents can be labeled 
as {mi=f, m2=l 7 rri3=T, rri4=f}. Each of the .4.(12) and A(34) atoms has an atomic spin Sf 2 ,S£ 4 = or 1. As a 
consequence, the s-wave molecular state has a total molecular spin S££ = 0, 1, or 2. By the definition of the 9-j 
symbols [4l| , the probability amplitude for the occurrence of fermion-antifermion spin states of S{{ and S l J 2 in a state 
with atom spins S^ 2 and and molecular spin is 



/ off ell . qAAi qA qA , qAA\ 
Wl4 D 23> '-'tot \ D 12 D 43^ d U 



J tot 



si s 2 5f 2 



(2S[{ + l)(2Sg + l)(2Sf 2 + l)(2S£ 3 + l){ s 4 _ *3 S&_ } , (60) 

off ell qAA 
.°14 °23 "-'tot . 

where si=s 2 =s 3 =s 4 =l/2. The probability \(S(( S l J 3 ; S^flS^S^; S^f}\ 2 for different atomic spins of S^ 2 and 
combining into different S{{ and S l J 3 , for a fixed total molecular spin Stot, are given in Table III. 
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TABLE III. Probabilities for different atomic spins of S12 and S34, combining into S{{ and S23, for a fixed total molecular 
spin S&t- 





qA 
^12 


qA 


qII 


oil 
'-'23 


probability 


QAA 
°tot 


= 2 


1 


1 


1 


1 


1 


QAA 
°tot 


= 1 


1 


1 


1 





1/2 






1 


1 





1 


1/2 


QAA 
J tot 


= 1 


1 
1 
1 







1 
1 




1 



1 


1/2 
1/4 
1/4 


QAA 
°tot 


= 1 







1 
1 
1 


1 
1 




1 



1 


1/2 
1/4 
1/4 


°tot 


= 


1 
1 


1 
1 


1 




1 




1/4 
3/4 


QAA 
^tot 


= 










1 




1 




3/4 
1/4 



TABLE IV. The s-wave (/x + e )-(e + /x ) molecular state energies and annihilation lifetimes in different Stcit > S12, an d 534 spin 
configurations. 



Molecular state 


Energy(o^) 


QAA 
°tot 


qA 


QA 
°34 


Annihilation 
Lifetime r an n(sec) 


Is 


-33.5 


2 


1 


1 


0.861 xl0" s 






1 


1 


1 


0.154xl0~ 10 






1 


1 





0.308 xl0~ 10 






1 





1 


0.308 xlO -10 









1 


1 


0.103xl0~ 10 















0.308 xlO" 10 


2s 


-9.30 


2 
1 
1 
1 




1 
1 
1 



1 




1 
1 



1 
1 




0.405x10"'' 
0.726xl0" 10 
0.145xl0~ 9 
0.145xl0~ 9 
0.484xl0~ 10 
0.145xl0" 9 



By taking into account spin probabilities^ the total rate of annihilation for matter-antimatter molecules consisting 
of 4 fermions with total molecular spin S££, initial atomic spins 5*^ and S^ 4 is 

w^S&S&SA*) = ]T \{SuS 23 ; Sj£)\S&S&; S^flw^ + w lT s J, (61) 

Sl4<S23 

where for lepton annihilations w 1 ^ is given by Eqs. (|58[) and f)59[) . The annihilation lifetime is then given by 



^S&S-uSfo?) - (qA qA qAA , (62) 



The rate of annihilation of the (/i + e )-(fi + e + ) molecule in its ns states can be calculated from Eq. (jfTTj) by identifying 
/=/i and l=e. The results for the annihilation lifetimes are shown in Table IV where we list the (/i + e~)-(e + /j~) 
molecular state energies, the molecular spins S^ti * ne atomic spins S^ 2 , £34 and the annihilation lifetimes T aim - In 
calculating the annihilation lifetimes in their physical units, we have used Table I to convert atomic units to physical 
units. 

The (/z + e _ )-(e + /i _ ) molecular states with spin S££=2 have the longest annihilation lifetimes, corresponding to 
lepton-antilepton pairs in their spin triplet states. The molecular Is and 2s S^f = 2 states have annihilation lifetimes 
of 0.861 x 10~ 8 and 0.405 x 10 -7 sec, respectively. Their relatively long annihilation lifetimes may make them accessible 
for experimental observations. 
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C. Annihilation of Hadron Constituents 



Hadrons are composite particles consisting of quarks and/or antiquarks whose quantized masses are governed 
by non-perturbative quantum chromodynamics. As a consequence, there are no simple selection rules for hadron 
annihilation similar to those for the annihilation of leptons. 

The lightest quantized hadrons are pions. Because of the difference in their masses, the annihilation of heavier 
hadron pairs such as pp and K + K~ differ from the annihilation of 7r + 7r _ . The annihilations of pp and K + K~ pairs 
proceed through strong interactions as many pairs of pions and other hadrons can be produced. In contrast, a tt + tt~ 
pair can annihilate through strong interactions only when its center-of-mass energy y/s exceeds the threshold of Am^, 
with the production of an additional pion pair. For pions in a matter-antimatter molecule, the momentum of the pion 
is of order am, and it has energy much below the strong-interaction annihilation threshold. We can infer that pions 
in a matter-antimatter molecule annihilate predominantly through the electromagnetic interaction with the emission 
of photons or dileptons. 




FIG. 11. pp total cross section in mb as a function of pi a b 



We shall first discuss the annihilation of heavy hadrons such as pp and K + K~ . We envisage that these hadrons an- 
nihilate essentially through a geometrical consideration depicting the occurrence of annihilations within a geometrical 
area jtq, with important initial state interactions that lead to a Gamow-factor 1/v type enhancement at low energies 
[42l l43f . We can therefore parametrized the annihilation cross section as 

0"ann = — , (63) 
V 

where v is the relative velocity between the colliding hadrons. From the PDG data J44|, the total and elastic pp cross 
sections in Fig. [TT] obey the following relationship: 



where 

Pub 



(65) 



As pp inelastic cross section is the same as the pp annihilation cross section, the experimental data gives ero = 50 
mb, and 

S"ann« = °0 = 50 mb. (66) 

For other systems, we can use the additive quark model p5l . |4(| to infer that 

ait « < x n\, (67) 
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where and are the number of quarks in h and h. We get 



^ ^50mb. (68) 



Consequently, the rate of annihilation of constituents 1 and 4 through strong interactions per unit time is 



w 



hh _ n q' l q 



14 " g 



50mbxPi4, {hh} = {pp} ot {K + K~}. (69) 



In the (pe )-{e + p) molecule, there are two leptons and two baryons. The spins of the four fermions are coupled 
together and we have 

w aa n(S&s*s£*) = u&+ \(ffls£*- ] stf\s&s&s£)\*wij + . (70) 

qPP qe+e — 

In the [K + e~)-(e Jr K~) molecule, there is a lepton pair and a K-K pair. As the kaons have spin zero, the molecular 
spin Stat comes only from the leptons and we have 

U'ann(Stot) = ™U "~ + wf^~ . (71) 

In practice, for the molecular ns states we have considered, the rate of hadron annihilation is much greater than the 
rate of lepton annihilation, if the hadrons can annihilate through strong interactions. In this case, because of the 
dominance of the hadron annihilation through strong interactions, the total annihilation rate is essentially independent 
of the spin of the molecule Stot- 

For the case involving 7r + and 7r _ constituents, annihilation cannot proceed through strong interactions as the 
pion energies are below the inelastic threshold. The pion pair can annihilate into photons and dileptons. Because 
the tt + -it~ system has the same total angular momentum and parity quantum numbers as the spin singlet state of a 
lepton pair, and the Feynman diagrams for the emission of two photons in QED for I + I — > 2j and 7r + + n~ — > 27 
have the same structure, we can approximate the 7r + +7r~— > 27 cross section to be the same form as the spin-singlet 
I + I — > 27 cross section in Eq. (|56p , 

~ t (-^) ■ (72) 

The cross section for pion annihilation into dileptons [as given by Eq. (14.33) of Ref. [45| ] is of order a 3 for pions 
with momentum p ~ am w . We can neglect the dilepton contribution from 7r + +7r~— > l+l in the present estimate. 
Consequently, the rate of annihilating pion constituents 1 and 4 by electromagnetic interactions per unit time is 

<4 +7r " = 7T (-^) X P U . (73) 

In the (7r + e _ )-(e + 7r _ ) molecule, there is a lepton pair and a tt + -t:~ pair. As the pions have spin zero, the molecular 
spin comes only from the leptons. For the (7r + e~)-(e + 7r~) molecule, we have 

^ann(Stot)=^4 +,r_ +^lf ■ (74) 

Because the annihilation of both the lepton pair and the pion pairs are electromagnetic in origin, they are comparable 
in magnitude. The annihilation rate depends on the spin of the molecule Stot- 

In Table V, we list the 1=0 molecular states, their energies, and their lifetime for (pe~)-(e + p), (p^~)-{^ + p), 
{K + e _ )-(e + A' _ ), (A' + /i _ )-(/i + A' _ ), and (7r + e~)-(e + 7r~). We observe that the annihilation lifetimes are short for the 
(pe + )-(e + p) states, of the order of 10~ 15 -10~ 18 sec, increasing to order 10~ 10 -10 -11 sec for (7r + e~)-(e + 7r~) states. 



XI. SUMMARY AND DISCUSSIONS 



To examine the stability of matter-antimatter molecules with constituents (mf rh\ ) and mi>ra2, we reduce 

the four-body problem into a simpler two-body problem. This is achieved by breaking up the four-body Hamiltonian 
into the unperturbed Hamiltonians of two atoms A(m^m^) and A(m^~mj~), plus residual interactions and the kinetic 
energies of the atoms. The unperturbed Hamiltonians of the two atoms can be solved exactly. Molecular states can 
be constructed by using the atoms A and A as simple building blocks. The interaction potential V(r) between the 
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TABLE V. Molecular Is state energies and annihilation lifetimes for (pe )-(e + p), (pfi )-(/i + p), (K + e )-(e + K ), (K + fi )- 
{H + K~) and (7r + e")-(e + 7r"). 



System 


State 


Enereyfa^u) 




Annihilation 












Lifetime(sec) 




(ve~)-(e~^v) 


Is 


-306.208 




0.549xl0" iY 






2s 


-92.565 




0.347 xlO" 17 






3s 


-43.255 




0.107xl0" 15 






4s 


-24.595 




0.244xl0" 15 






5s 


-15.595 




0.464xl0~ 15 






6s 


-4.094 




0.259X10 -14 




(p/i _ )-(/i + p) 


Is 


-0.540 


0.532xl0" iY 


(K+e-)-(e + K-) 


Is 


-160.601 




0.853xlO" Ib 






2s 


-48.108 




0.526xl0~ 15 






3s 


-22.141 




0.166xl0" 14 






4s 


-12.316 




0.372xl0~ 14 






5s 


-7.581 




0.715X10 -14 






6s 


-4.947 




0.122xl0" 13 




{K + ^')-{^K') 


Is 


-0.0239 


0.154xlO~ Ib 


(7r + e-)-(e+7r-) 


Is 


-44.625 


(Stot=0) 0.592x10" 


ii 








(S tot =l) 0.622x10" 


ii 




2s 


-12.698 


(Stot=0) 0.298x10" 


10 








(S to t=l) 0.392x10" 


10 




3s 


-5.341 


(S tot =0) 0.616x10" 


10 








(S'tot=l) 0.121x10" 
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atoms is then the sum of the direct potential Vdi r (?") arising from the interaction between the constituents and the 
polarization potential V po \(r) arising from the virtual excitation of the atomic states in the presence of the other 
atom. The eigenenergies of the molecular states can be obtained by quantizing the four-particle Hamiltonian (22j . 

The Coulomb four-body system in atomic units of A and A atoms is completely characterized by mi/m2. Con- 
sequently, two different four-body systems with the same mx/m2 will have the same molecular state eigenenergies 
and eigenfunctions in atomic units. We can infer the stability of a molecule by studying the change of the molecular 
eigenenergies as a function of toi/to2- 

The effective reduced mass of A(rrii m^)-A(m^"mi) for molecular motion is large when mi/ni2 3> 1 and decreases as 
mi/ rri2 approaches unity. The relative importance of the direct and polarization potential also changes with m\/m2- 
For a matter- antimatter molecule with m\Jm,2 3> 1, we find that the direct potential dominates in regions of r < a 
and gives rise to deeply bound molecular states. As m\/m2 approaches unity, the magnitude of Vdir decreases and 
the polarization potential by itself is too weak to hold a bound state. As a consequence, the state energies (in atomic 
units) rises and comes up to the continuum as mi/m2 approaches the unity limit. 

We find that matter-antimatter molecules possess bound states if mi/m2 is greater than about 4. This stability 
condition suggests that the binding of matter-antimatter molecules is a rather common phenomenon. This molecular 
stability condition is satisfied, and many bound molecular states of different quantum numbers are found, in many four- 
body systems: (/i + e~)-(e+/i~), (7r + e~)-(e + 7r~), (A' + e~)-(e + A'~), (pe~)-(e + p), (p^i~)-(fi + p), and (K + [j,~)-(n + K~). 
Bound molecular states are not found in (r + p)-(pT _ ) which has m\jm2 = 1.89. 

When one applies the stability condition to the (e + e~)-(e~e + ) system, one may naively infer at first that the 
(e + e - )-(e~e + ) system will not hold a bound state. On the other hand, bound (e + e~) 2 molecular state has been 
experimentally observed [llj], and the binding energy has been calculated theoretically to be 0.016 a.u. relative to two 
(e + e~) atoms Q. It should however be realized that the stability condition we have obtained applies to four-body 
systems with distinguishable constituents without identical particles. For systems with identical particles such as the 
(e + e~) 2 system, it is necessary to take into account the antisymmetry of the many body wave function with respect 
to the exchange of the pair of identical particles. Depending on the spin symmetry of the identical particle pair in 
question, the symmetry or antisymmetry with respect to the spatial exchange of the pair will lead to a lowering or 
raising of the energy of the state. The identical particles in the bound (e + e~) 2 system leading to the bound state 
have been selected to be spatially symmetric states [l|, , which corresponds to spin-antisymmetric with respect to 
the exchange of the pair of identical particles Q- As a consequence, their spatial symmetry lowers the state energy 
relative to the state energy when there is no such a symmetry. The small value of the theoretical binding energy 
(0.016 a.u. 0) suggests the occurrence of such a lowering of the energy from the unbound to the bound energy region. 
In order to confirm this suggestion, it will be of interest in future work to extend the present formulation to include 
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the case with identical particles and their exchange symmetries. How the additional symmetry considerations may 
modify the stability of those molecules with identical particles is worthy of future investigations. The states we have 
obtained may not contain all the correlations. Additional correlations may be included by adding correlation factors 
and using the solution of the present investigation as doorway states. Future addition of correlations superimposing 
on the states we have obtained will be of interest. 

We can divide the annihilation lifetimes of matter-antimatter molecules into different groups that correlate with 
the types of constituents. Those molecules constructed from different leptons have the longest annihilation lifetimes, 
of the order of 10 _8 -10 -11 sec, depending on the spins of the molecular state. The second group involves leptons 
and pions in which the pions cannot annihilate through strong interactions and can annihilate only through the 
electromagnetic interactions. The annihilation lifetimes are of order 10 -10 -10 -11 sec. Molecular states containing 
kaons have annihilation lifetimes of order 10~ 14 -10 -16 sec while those containing proton and antiproton 10 _15 -10 -18 
sec. The relatively long annihilation lifetimes for leptonic (fi + e~)-(e + fj,~) molecules may make them accessible for 
experimental detection. 

We have examined only molecular states constructed from the AA family in which the building-block atoms Aq and 
Aq are in their ground states. We can likewise construct in future work states in the AA family in which A and A are 
in their excited states. These states will have different interatomic interactions, obey different stability conditions, 
and have different properties with regard to annihilation and production. In another future direction of extension, we 
can also construct molecular states based on the MM family, using states of M(mffh^) and M(m^m^) as building 
blocks. The AA molecular states and the MM molecular states have different topological structures and properties. 
Molecular states built on different branches of the family tree will likely retain some of their family characteristics. 
These possibilities bring into focus the rich degrees of freedom and the vast varieties of states that need to be sorted 
out in the Coulomb four-body problem associated with matter-antimatter molecules. 

In addition to the problem of molecular states as a structure problem investigated here, future theoretical and 
experimental studies should also be directed to the question of production and detection methods from reaction 
points of view. While the observation of new matter-antimatter molecules may be a difficult task, the prospect of 
classifying the (m^ fn^) system into the genealogy of families and basic building blocks, if it is at all possible, 
will bring us to a better understanding of the complexity of the spectrum that is associated with the complicated 
four-body problem opened up by the pioneering work of Wheeler 
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Appendix A: Evaluation of the Matrix Element 

{A x Ay\v(r jk )\A A } 

The evaluation of the excitation transition matrix element requires first the Fourier transform of the transition 
density. The transition density for the excitation from the ground 1 100) state to the excited \nlm) state is 



N 10 N nl f 2 
■ exp 



V47T 

Making the scale transformation r = r 1 J (n + 1) = r'/3, we have 



^}(^K^>w>. (ad 

2nao J V nao J nao 



V47T I 2na J \nao J \na 



PnimMr) = ^T^cxp pL*_ 2 Y* m {6^). (A2) 



We can expand L"(f3x) as a sum over Laguerre polynomials L°L k (.t) as given in Eq. (22.12.6), page 785 in [3t 

L?(J3x) = J2 ( j t a ) ^'""C 1 ~ PT L *-M- (A3) 



3 v > — — ~ — ™— — 

3 

K 

The transition density becomes 

n-2 



woo(r) = ^" 2 " K ( 1 - /^^WOoM, (A4) 



^ -=o 



where 

^■nim,ioo( r ) 



2r' 1 3 / 2/ \ 



m \ i 
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Equation (|A4[) is a sum of many terms, each of which is a hydrogen wave function. The Fourier transform of 

^ni™,ioo( r ) is s ivcn b y 

^nlm,10o(P) — J ^ r e ^ r ^nlm,10o( r ) 

= ^[Yrm (^)7 * 2 d XJl (^yy) ^xLU_ K (*) . (A5) 

Using the generating function of the Laguerre and Genegnbauer polynomials, the Fourier transform can be carried 
out |47[ , and we obtain 

4« („\ A-wmf^v* f na o\ 3 npa (n + l) 2 2 6 (n - k) 

^-^f i^Ti-^D - (A6) 



V(npa ) 2 + (n + l) 2 

Utilizing this result, we can write down the Fourier transform of the transition density in the form 

P„l m ,10o(P) = £nlm,100(p)K m (p)]*, (A7) 

where i?ni m ,ioo(? 5 ) is given by Eq. (J3SJ) . The transition matrix element becomes, 
(A x B x \Y, k \0O){T) = 8 J j^e^ r R^ 1Q0 (f A (jk)p)Y* n (p) 

x Rn' l-m, 100 (Jb {jk)p)Y?_ m (p)v jk (p) , 
V 2 dp ^ 



(27T) 
2?r 1 



(fAtik)p)R nll _ mA00 (f B (jk)p)v 3k (p) 



xjd^Jdn e* p - r Y* m (p)Yi*- m (p)- (A8) 

-1 

The angular integral can be carried out and we obtain 

A{P , r) = / *, y lTO (ff p ) l? m (* p ) = ? " ^ <** m ^ °>' (A9) 

J J \-\3o(pr) +J2(pr)\ (form = l), 

o-i v 

where jz(x) are spherical Bessel function. This gives Eq. (|36|) in Section VLA. 
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